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The staid subject of exact static spherically symmetric perfect fluid solutions of Einstein's equa- 
tions has been reinvigorated in the last decade. We now have several solution generating techniques 
which give rise to new exact solutions. Here the Einstein static Universe is transformed into a 
physically acceptable solution the properties of which are examined in detail. The emphasis here is 
on the importance of the integration constants that these generating techniques introduce. 
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I. INTRODUCTION 

A decade ago, as a demonstration of the use of com- 
puter algebra [l| , it was shown that few of the alleged ex- 
act static spherically symmetric perfect fluid solutions of 
Einstein's equations were in fact correct and even fewer 
made physical sense Since then this field of study 
has been reinvigorated with the development of several 
solution generating techniques which give rise to new ex- 
act solutions. These techniques are based on the pressure 
isotropy condition, either looked at as a differential equa- 
tion, following Wyman [3i] , or by way of invariance prop- 
erties, following Buchdahl For example, following the 
work of Rahman and Visser Q , I showed that several of 
the acceptable solutions given in [l| follow from a simple 
algorithm \d\ which in fact generates an infinite number of 
physically acceptable solutions. Since then several works 
of interest have appeared: 0, B, 0, O, [Hi- [H and 
(Tst . It is remarkable that such a staid old subject has 
bounced back to life. 

Here I do not add to these generating techniques but 
rather make use of one of them to do what they are in- 
tended to do: generate a new physically interesting exact 
solution of Einstein's equations. The example provided 
serves to emphasize the importance of the integration 
constant that the generating technique introduces. 



II. GENERATING TECHNIQUE 

Every spacetimc Ai with metric [ij] 
dr^ 



where 
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2m(r) 



+ rW-e2(*W+«W)dt2, (1) 



where dfl^ = dd^ + s\n{9Yd(jP , is an exact perfect fluid 
solution of Einstein's equations as long as 



and 



m — T. 

gj tt(r)(ir 
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= and the C„ are constants. The functions a(r), b(r) 
and c(r) are given by 
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The algorithm can be executed subject to the specifica- 
tion of the function 4> (as well as smoothness and bound- 
ary conditions f6l|) and the constants C„ [Tsj . 

The procedure I consider here will assume C2 = and 
so C3 is disposable (it can be absorbed into the scale of 
t). Call these spacetimes N. Further, I will write the 
spacetimes in the form 



dsQ 



dr^ 



2M(r) 



(8) 



where M, as with m in A^, is constructed so as to make 
O a perfect fluid. That is, for M , ([5]) is replaced by 



2r2($" + $' 2) _ 3(j.$- + 1) 4r2(x" - X^) - 6rx' 



"(r^' + 1-1- 2rx 



(9) 



and ([6]) by 

- _ r(r($" + ^' ^) - + 2{rx - rx ^ - x )) 
r$' -I- l-f 2rx 

in Unlike m however, M is no longer the effective 
gravitational mass @. 

It is important to note that O is not a conformal trans- 
formation of A/" (due to the restrictions on M) and it is no 
more general than J\f as it is merely a coordinate trans- 
formation of A/" [1^ . We refer to the case x = as the 
"seed" of the spacetimes O. The usefulness of the form 
([5]) derives from the fact that we can clearly recognize 
the seed. 
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III. $ = 

The simplest seed for O is $ = C where C is a constant 
which, by choice of scale for t, we can set to zero. It 
foUows that m = Cir^ and the seed is simply the Einstein 
static Universe. (A cosmological constant A = 2Ci can 
be introduced to give zero pressure, but this is not done 
here |l3)- Given this seed, the regularity conditions on 
X are [ij 



|X(0)|<(X3, x(0) = 



(11) 



and so the simplest non-trivial form of x satisfying these 
conditions is 



X = C4 + Car 



(12) 



where C4 and C5 7^ are constants [l8j. Since the con- 
stant C4 simply scales the energy density and pressure by 
1/e^'''', without loss in physical generality we set C4 = 0. 
Since C5 can be absorbed into a redefinition of r (and a 
rescaling of the as yet to be chosen constant Ci) we set 
C5 = 1 so that without any loss in physical generality we 
take X = r~. We now have 



(13) 



where R = vr+4r^ and erf is the error function. 
Whereas we could of course continue our discussion in 
terms of the coordinates used in ([8]), we now revert to 
more traditional coordinates. 

Under the coordinate transformation r = r we now 
have 



2^ 



-de 



where 



with 



{J + 2r2(V2^ £ - Ci)) (1 + nf 



£ = erf( 



and 



n = W(2r2) 



(14) 
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where W is the Lambert W function [l9|. As shown in 
FIG. [U the constant Ci plays the central role regarding 
the physical acceptability of these solutions. For 



(19) 



p vanishes while p > which is physically unacceptable. 
For 



(20) 



the solution is global, not isolated, and p and p ^ only 
as r — > c». Finally, for 

1 



V27re <Ci < 2 + V27re erf(^) (21) 

v2 

the pressure p vanishes at finite r > and the solutions 
match onto a vacuum exterior by way of continuity of 
the effective gravitational mass. The density contrast at 
the boundary increases as Ci increases. Throughout the 




FIG. 1: The abscissa is r and the ordinate is Ci. Top down 
the curves show p =0, p = 0, p = and p = Q. The curves 
p = and p = intersect the ordinate at 2 + VSvre erf (-^) 

and \/27re erf(-^) respectively. They are both asymptotic 
to Ci = ^27re. 



physically acceptable distributions p and p are monotone 
decreasing and the adiabatic sound speed Vs is sublumi- 
nal. Some properties are shown in FIG. [21 

We record here in explicit form the essential physi- 
cal elements of the solutions: The effective gravitational 
mass is given by 



M = 



(22) 



with T given by (jlSp . the energy density is given by 
2/Cr2 (Ci - V2^ £) - 3 £ 



Stt p 



(23) 
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FIG. 2: Some internal properties of the solutions. Top down 
at the origin the curves give p/2, p, Vs and 2M/r where Vs is 
the adiabatic sound speed and M the effective gravitational 
mass. The main image is for asymptotic case Ci = V2Tve and 
the insert is for Ci = 4.24 which terminates at p = where 
r ~ 0.7914. 



(28) 



Finally, the square of the adiabatic sound speed is given 

by 



where 



and 



V = J' 



,2nP (^/2^£ -Ci) +S 
2Tt^ (Ci - £)+U 



n={j' ^i){J -1){J + l)J, 



T={J^-J' + 6)(J -1){J+ 1) J, 



U 



- Xhj"^ + 30. 



(29) 
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IV. DISCUSSION 



where 



and 



and the isotropic pressure is given by 

2-Pr2 {^/2^E-Ci) + Q 



where 



V = {ij^ - 1){J^ + 1) 



(24) 



(25) 



(26) 
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The Einstein static Universe has been transformed 
into a class of physically acceptable static fluid spheres 
whose physical properties have been written out in ex- 
plicit form. The technique is but a coordinate transfor- 
mation of one discussed previously @] but it allows a 
clear understanding how various spacetimes can be in- 
terrelated [201. 



Acknowledgments 

This work was supported by a grant from the Natural 
Sciences and Engineering Research Council of Canada. 
Portions of this work were made possible by use of 
GRTensorll [H. 



[*] Electronic Address: lake@astro.queensu.ca 
[1] See M. S. R. Delgaty and K. Lake, Computer Physics 
Communications 115, 395 (1998) arXiv:gr-qc/9809013 . 
A partial list of corrections can be found at 
http : //grtensor .phy . queensu. ca/ solutions/ 
[2] The conditions used in [l|] were: (i) isotropy of the pres- 
sure (otherwise any metric is a "solution"), (ii) regular- 
ity of associated invariants at the origin, (iii) positivity 
of the pressure and energy density, (iv) vanishing of the 
pressure at a finite boundary, (v) monotone decrease of 



the pressure and energy density to the boundary and (vi) 
subluminal adiabatic sound speed. In addition to these, 
a monotone decrease in the subluminal adiabatic sound 
speed was considered desirable. Probably only issues re- 
garding the "sound speed" would offer room for debate. 

[3] M. Wyman, Phys. Rev 75, 1930 (1949). 

[4] H. A. Buchdahl, Phys. Rev. 116, 1027 (1959). 

[5] S. Rahman and M. Visser, Class. Quant. Grav. 19, 935 
(2002) arXiv:gr-qc/0103065 . 

[6] K. Lake, Phys. Rev. D 67 (2003) 104015 



4 



[arXiv:gr-qc/0209104'. 
[7] D. Martin and M. Visser, Phys. Rev. D 69 (2004) 104028 

[arXiv:gr-qc/0306109; . 
[8] P. Boonserm, M. Visser, and S. Weinfurtner, Phys. Rev. 

D 71 (2005) 124037 arXiv:gr-qc/0503007 . 
[9] P. Boonserm, M. Visser, and S. Weinfurtner, Phys. Rev. 
D 76 (2007) 044024 arXiv:gr-qc/0607001 . 
[10] P. Boonserm, "Some exact solutions in general relativ- 
ity" , MSc thesis, Victoria University of Wellington, 2005 
[arXiv:gr-qc/0610149] . 
[11] P. Boonserm, M. Visser and S. Weinfurtner, Jour- 
nal of Physics: Conference Series 68 (2007) 012055, 
arXiv:gr-qc/0609088;. 
[12] P. Boonserm, M. Visser and S. Weinfurtner, "Solu- 
tion generating theorems: Perfect fluid spheres and the 
TOV equation", |arXiv:gr-qc/0609099 (Marcel Gross- 
mann 11). 

[13] P. Boonserm and M. Visser, Int. J. Mod. Phys. D 17, 

135 (2008), arXiv:0707.0146vl]. 
[14] I use geometrical units and usually designate functional 

dependence only on the first appearance of a function. 

Throughout p desinates the energy density and p the 

isotropic pressure. 
[15] Up to notation and a coordinate transformation the al- 



gorithm given is equivalent to Theorems 1 and 2 in p^ . 

[16] Define re^'^^^ = r and <l>(r) -I- x(t') = V'('') ^'id write 
^j-j. = 1/(1 — 2m(r)/r). Then m(r) is given by (0 with <1> 
replaced by tp and r by r. 

[17] For a recent discussion of perfect fiuid spheres with A see 
C. Bohmer and G. Fodor, Phys. Rev. D 77 (2008) 064008 
[arXiv:0711.1450 . 

[18] This solution is distinct from but in a sense compli- 
mentary to the solution generated by y = and <1> = 
Ci -f Csr^; Kuch2 111 in the notation of Q[. 

[19] W is defined by the condition W(a:)e^'^) ^ x so that 
gtt in p4p can be written in the equivalent form —e^. 
For a discussion of W see, for example, R.M. Corless, 
G.H. Gonnet, D.E.G. Hare, D.J. Jeffrey, and D.E. Knuth, 
Advances in Computational Mathematics 5, 329 (1996). 

[20] Closed form exact solutions have also been found for 
other choices of "1> and x by Alex Klotz, Cedric Grenon 
and Pascal Elahi. 

[21] This is a package which runs within Maple. It is entirely 
distinct from packages distributed with Maple and must 
be obtained independently. The GRTensorll software and 
documentation is distributed freely on the World-Wide- 
Web from the address htt p:/ /grtensor. org\ 



